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ABSTRACT 


The United States Air Force uses a nonlinear programming model to assess the 
utilization of weapons and sorties needed to achieve a maximum value of destroved 
targets in a multi-period, Theater-Level conflict. The current model is modified bv con- 
straining the consumption of weapons. Alternate objective functions are introduced. 
Their meaning and influence on the optimization 1s compared. An increase in the worth 
of destroyed targets is gained if the model can more flexibly utilize weapons than 1s 
currently the case. The optimization can be further improved if all time periods are 
considered simultaneously while assigning sorties to targets, rather than the current 


myopic approach. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While everv effort has been made, within the 
ume available, to ensure that the programs are free of computational and logic errors 
thev cannot be considered validated. Any application of these programs without addi- 


tional verification is at the risk of the user. 
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l. INTRODUCTION 


In 1988 the United States Air Force purchased over $ 2 billion worth of weapons for 
use in different theaters around the world. The projected need for the quantity of differ- 
ent Weapon types is based on an annual Nonnuclear-Weapon Consumables Analvsis 
(NCAA) performed by the Directorate of Plans, USAF [Ref. 1]. Unlike other services, 
the USAF relies widely on mathematical programming models in order to optimize the 
allocation of weapons. 

In 1974 RAND developed a nonlinear programming model that optimizes the 

number of different sortie types assigned to several target tvpes bv maximizing the mill- 
tarv worth of killed targets [Ref. 2: p. 5]. Since each target tvpe was given a different 
target value, the model attempts to assign sorties to maximum value targets first. To 
avold an undesired concentration of sortie allocations to a few or even one target tvpe. 
a nonlinear objective function was introduced. Within the model onlv the number of 
available targets and sorties are constrained. The expenditure of weapons is not consid- 
E Ihe number of targets one Sortie is able to destroy is expressed by an effectiveness 
parameter that depends onlv on sortie and target tvpe. 
The required input data structure for the RAND-model is a simplification of the 
much more complex data base contained in the Joint Munitions Effectiveness Manual 
(JMEM) used bv USAF. The JMEM data base determines effectiveness as a function 
of weather and mission profile (tactic) as well as type of aircraft and type of target. In 
the current operation a model called SELECTOR sorts the JMEM data base so that for 
each sortie-target type combination, all feasible tactics are ordered from the most to the 
least cost-effective, including the cost of aircraft attrition. This list 1s referred to as the 
Preferred Weapon List. 

The data in the Preferred Weapon List must be reduced to input parameters de- 
pending only on sortie and target type as mentioned earlier. This 1s basically done by 
selecting the most cost-effective tactic from the list feasible for weather situations con- 
sidered in the model. After the optimization has determined the optimal number of 
sorties assigned to different targets, the number of remaining targets and the expenditure 
of weapons is evaluated. This process is repeated in subsequent time periods with a new 
inventory of sorties and also by recording the remaining number of active targets and 


weapons available. In this way, tactical changes in a given scenario over time are 





considered bv optimizing sequentially for discrete time periods. This process is 
accomplished in one programming model and is called the HEAVY ATTACK model. 
The USAF interest 1s mainlv in the consumption of weapons utilized over all time 
periods. 

The objectives of this Thesis are to include a weapons constraint in the 
RAND-model and to investigate alternatives to the currently used objective function. 
In addition, the RAND-model is expanded so that more available information is in- 
cluded in the optinuzation in order to gain a higher total military worth of killed targets 
than is currently achieved. Therefore, the consumption of weapons used bv less cost- 
effective tactics 1s investigated when other weapons, used by the most cost-effective 
tactic, are exhausted. As a final consideration, one global optimization over all time pe- 
riods is compared to the current sequential optimization method. Global optimization 
achieves a higher overall worth of killed targets. However. gaining a higher militarv 
worth of killed targets serves onlv as an aid in analvzing the predicted need of weapons. 
The value of the revisions suggested in this Thesis have to be measured on their ability 


to satisfv the demands of the USAF and simultaneouslv meet budget constraints. 


П. BASIC STRUCTURE OF HEAVY ATTACK 


A. THE ORIGINAL RAND - NODEL 

In 1974 RAND developed a nonlinear programming model whose objective was to 
determine the optimal number of sorties of type i assigned to targets of type j by maxi- 
mizing the total military value of destroyed targets. The relationship between an as- 
signed sortie and a target kill is established by introducing "sortie effectiveness” D МЕ 
parameter Ё, , defines the average number of kills that one sortie of type 1 will achieve 


when it is assigned to targets of tvpe j. 


Definition of index 


1 sortie tvpe 
j target tvpe 
Parameter 


T, total number of tvpe j targets available at the beginning of a time period 
V, military worth of tvpe j target 
S 


t 


total number of tvpe 1 sorties available 


E average number of type j targets killed by one type 1 sortie 
Variables 


Б 


jj number of type i sorties assigned to type j targets 


Моде! 


Max z= > 1; x дув, х sx) 
j i 
5 
Эр R vi 
7 


ДУ) < Т, у) 
і 

2,53, € cx ? SA, vi 

jeJ J 

where J i$ a subset of all targets OT DEI OO eT 

О = 54; ; V 


(220. S.V) is a concave function that approaches | for large arguments. [he 
RAND - model (and HEAVY ATTACK) utilizes a specific analvtic from that will be 
examined in detail later. The recipe constraints ) SX, = ¢ x >5Л limit the numer 
of sorties of tvpe 1 which are assigned to a list of targets bva fraction of the total number 
of sorties of tvpe 1. Since these constraints are not used bv the USAF in their current 


weapon analvsis, this inequality will omitted from now on in the Thesis. 


B. THE ROLE OF SELECTOR 
Based on the information contained in the JMEM the effectiveness of a sortie de- 


pends on sortie type, target type, Weapon type , weather and tactics or mission profile. 


Definition of index 


1 sortie type 

] target type 

k Weapon type 

ту weatherband index 
Г index for used tactic 


Definition of parameter 


Е,,,, number of tvpe | targets killed by one type 1 sortie using tactic r 
in weatherband w 

B,,,, number of weapons carried bv one type 1 sortie which is assigned to tvpe j 
target in weatherband w and using tactic r 

K 4 ipe of weapon which is loaded on sortie 1 and will be deployed to target j 


bv using tactic r in weatherband w 


The JMEM data have too many subscripts to match the required input data structure 
of the RAND - model. The number of subscripts of a sortie needs to be reduced so that 
E, depends only on sortie and target type. The first part of the task of reducing the 
number of subscripts from 4 to 2 is accomplished by the sorting program SELECTOR. 
The output data of SELECTOR - referred to as Preferred Weapon List - contains for 


each different sortie - target tvpe combination five distinct items: 


l. The worst weatherband in which a tactic can be used. 


لم 


The tvpes of weapons that can be allocated. 


оз 


Niietreicriwe COst-elliciency ol @ tactic given by Is Order on the list. 


4. The number of targets which can be killed by one sortie. 


Ca 


The number of weapons that can be carried by one sortie for each weapon type 
(mixes of Weapons are not considered). 


The data structure of the Preferred Weapon List, which will be used later for the aggre- 
gation of the input data E,, for the RAND - model, is illustrated by the following ex- 


ample: 


Subset of data from Preferred Weapon List 


1 ) I w I Е В, ЯЛЖ) 
Рр 22 1 4 3 0.137 4 
12 2 B 1 0.664 6 
П, B 2 17 1.580 2 
ША 9 4 3 17 1.000 2 


For example, the most cost-efficient and feasible tactic for weatherband w- 3 is 
tactic r=2. Tactic r=1 1s more cost-efficient because it 1s first on the Dist, but is only 
feasible in weatherband w=4 or higher. Weatherband w= 1 expresses best weather while 
weatherband w 6 represents the worst weather. Tactic r= 3 1s feasible (a tactic feasible 
in w is alwavs feasible in better weatherbands) but less cost-efficient than tactic r= 2. 


The given data can be represented in the following way: 


Table ]. JE. VALUES: \umiber of tarcets of 1: 1ط برعم‎ ce 
1 using tactic r in Weatherband w. 


Tit [wt | we? | wed | wea | wed [wee 
ojs e Ге oe oer oe oss 
ооо рою 150 зво зво | зво 
ооо Го | ъс Lus > 







Чаше 2. Вы, - WALUES: Number of weapons that are loaded on one sortie of 
tvpe 1 which is assigned to target tvpe } and using tactic r in weatherband 





Wable 3. K,,24-- VALUES: Type of weapon that is allocated to a sortie of tvpe 1 
which is assigned to a target of tvpe j and using tactic rin Weatherband w 





Since HEAVY ATTACK only considers the tactic at the top of the lst for each 
weatherband, and since weapon tvpe is implied by tactics. SELECTOR essentially re- 


duces the number of subscripts from 4 to 3. 


С. DETERMINATION OF E,, INHEAVY ATTACK 

An important assumption for HEAVY ATTACK in order to understand the logic 
behind the aggregation of E,, is that the weather is nor known at the time when sorties 
are assigned to targets. This leads to the condition that the effectiveness of a sortie and 
the consumption of weapons in a particular weatherband has to be proportional to the 
probability that this weather will occur. 

This probability 15 E in HEAVY ATTACK by a given distribution of 6 


distinct weatherbands: 


PR, = probability that weatherband w will occur at a certain time in the future, 


pec po». 6. 


Throughout this Thesis the following distribution 1s used: 


Table 4. WEATHER DISTRIBUTION IN HEAVY ATTACK: Probability that 
weatherband w occurs when sorties are allocated to targets. 


[wet | w=? | was | wea | wes | wae 






Since weatherband w=1 will never occur, the effectiveness for any sortie in this 
weatherband 1s irrelevant. Itis assumed that any weapon which is feasible for a certain 
sortie - target combination can be used in the weatherband determined by SELECTOR 
or in any better weather (higher weatherband). 

HEAVY ATTACK uses for each weatherband onlv the top weapon on Preferred 
Weapon List. This means that the model will allocate the most cost-efficient weapon 


feasible in each weatherband. Therefore the data set E. 


Г г.и 


can be reduced by the sub- 


script r such that: 


* 


. = the effectiveness of the most cost-efficient tactic in weatherband w. 


Table 5. EFFECTIVENESS OF THE MOST COST - EFFICIENT TACTIC: In 
each weatherband w the first effectiveness value in Table 1 greater than 
Zeros selected. 






Applving the same reasoning on the data set D,,,,, and K,,,,, Vields : 


ل 


* 


B;;y = number of weapons used by the most cost-efficient tactic in weatherband w, 


* 


Kia] = type of weapon used by the most cost-efficient tactic in weatherband w. 


Table 6. WEAPON LOAD OF THE MOST COST - EFFICIENT TACTIC: In 
each weatherband the first weapon load value in Table 2 greater than zero 
IS selected. 





Table 7. WEAPON TYPE OF THE MOST COST - EFFICIENT TACTIC: In 
each weatherband w the first weapon tvpe in Table 3 not equal to zero 16 
selected. 





Since each weatherband will occur with the probability PR, the averaged 


eflectivness must be 


= > PR, x Ej. = 0.240 


In general the process of obtaining B is a little more complicated than described 
above because HEAVY ATTACK is permitted to use tactics lower than first order when 
first order weapon tvpes have been exhausted. This can happen because HEAVY AT- 
TACK is actuallv a model of protracted war. First order tactics are preferred because 
they represent the most cost-effective tactic. The war may last for several periods (4 in 


this Thesis), and it is possible that certain tactics may not be feasible in later periods on 


account of weapon exhaustion. Suppose for example, that weapon tvpe 3 has been ex- 
hausted in a previous time period and is therefore no longer available. The top weapon 
for weatherband w=4. 5 or 6 is now weapon type I The new effectiveness values 
E 


га r. Ww 


ае. 


Table 8. E,,4,,,- VALUES AFTER WEAPON K=3 IS EXHAUSTED: Number 
of targets of type j killed by one sortie of type 1 using tactic r in 
weatherband w that 1s applicable. 


СТ Г [wat | w=? | wes [wea [wes [wee 
C ps [a | ХА | ХА NA рха рха | ХА 
По 3-9 [9 роз | 6s | 0.661 | 966 — 
ENNENENENEKNEKNEKNEEINEKIM 
p» а Го Го оо ee Ten 







Lsing the most cost-efficient tactic in each weatherband w gives the following ef- 


fectiveness values F` : 


AM 


Table 9. EFFECTIVENESS OF THE NEXT FEASIBLE COST - EFFICIENT 
IACTIC: In each weatherband w the first applicable effectiveness value 
in Table $ greater than zero is selected. 


Сте ре е рез ез ес 
зо | о 











which results in the averaged effectiveness: 


E = DPR, x Ej jw = 0.682. 


Note that the effectiveness has increased on account of the lack of weapon type 
k= 3! The SELECTOR output is ordered according to cost-effectivness (not effective- 
ness), SO it is quite possible that tactics far down in the Preferred Weapon List may ac- 
tually be quite effective. These tactics typically have high associated attrition, but 
attrition is not considered in HEAVY ATTACK once SELECTOR has done its job. 

Bv considering the same logic, it can be observed that the fourth order tactic on the 


Preferred Weapon List with £,,,, = 1.600 will never be used. This is because the third 


10 


order tactic uses the same weapon (in this case weapon type k= 17) in at least the same 
worst weatherband as tactic r= 4. 


D. TIME IN HEAVY ATTACK 

Once the effectiveness values £,, are evaluated, the required input data is available 
in order to optimize the number of sorties assigned to the different target tvpes. For 
most cases all targets are not Killed when the optimization is finished because of the 
constrained number of sorties in the RAND - model. As in a real war scenario, the 
outcome of a given attack will influence subsequent target consideration and planning. 
Only the targets that survived the previous attack will be reconsidered. Weapons are not 
resupplied and therefore may become exhausted. The current version of HEAVY AT- 
TACK may actually allocate ore weapons in a given period than are available at the 
beginning of the period. This is because there is no explicit constraint on Weapon usage. 
The deletion is currently done after each period by computing weapon usage after the 
optimization for the period is finished. However, a weapon will be deleted in the next 
period if it is exhausted at the end of the current period. 

There 1s no resupply of targets between periods in HEAVY ATTACK, although 
there is a facility for reconstituting targets that have already been killed. This will be 
discussed later. Aircraft are also not resupplied or even directly represented in HEAVY 
ATTACK: the number of sorties available during each period is a direct input. Each 
time period represents an attack which changes the input for the following time period. 

The fact that the importance of a target will change with time 1s represented in 
HEAVY ATTACK by the option of changing the military worth for each target type 
at the beginning of a new time period. Even though the military worth of a target 1s 
known in all future periods, the current sequential time optimization only “sces” the 
worth of a target for the current time period. Following from this “myopic” way of 
maximizing the military worth of killed targets it may happen that sorties are assigned 
in a time period to a target type when its military worth is relatively low. A "global (or 
overall) time optimization can be expected to achieve a higher military worth of killed 


targets. This is discussed later. 
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E. THE NONLINEAR MODEL IN HEAVY ATTACK 
The basic structure of the current model in HEAVY ATTACK for one time period 


number of type j targets available at the beginning of a time period 
number of dead tvpe | targets at the beginning of a time period 
military worth of type J] target during the current time period 
target - parameter for tvpe } target 

number of tvpe 1 sorties available for the current time period 
proportion of S, that can be assigned 


number of type i sorties that are assigned to type j targets 


number of type j targets killed in the current time period 


үсеп [Б 


Parameter 


PROP. 


Variables 


EL 


< 29 j 


РИ 


Model 
Mex 2 = > x KILL, 
J 


Sel. 


KILL, = 10 с, В), ЗЕ jX SA, 1 Y d 


where: 


z T; с, E 
До с. Р), 7 Бух sx) = 41 = 2) x ( ات‎ 2 т 


The above function is the same function as used by RAND [Ref. 2]. 


2 SXy € PROP, x 5, vi 
J 

"EU Dem г 27 
42 24 vij 


The nonlinear function /(T.. c, D, X E, x SX.) is of the same form as in the RAND 
- model. The number of targets of type j that are killed and the number of sorties of type 
1 are constrained. The consumption of weapons 1s not considered in the model itself. 
After the optimal numbers of sorties are determined by the optimization, the consump- 


tion of the different Weapon tvpes 1s evaluated by: 


( consumption of weapon ), = px x (228, x В. 
j w 


i 


ve 


where the sum 15 over allume such har m m 


F. TARGET RECONSTITUTION IN HEAVY ATTACK 

The abilitv to reconstitute killed targets is a common fact in a modern war. HEAVY 
ATTACK records the number and tvpe of targets as well as the time period when they 
are destroyed. After each optimization, it determines if targets can be reconstituted and 
evaluates the maximal number that are possible. A major task in this Thesis has been 
to determine the conditions under which reconstitution is allowed to happen by analyz- 
ing the responsible part of the HEAVY ATTACK source code. HEAVY ATTACK’s 


logic seems to be as outlined below: 


Definition of index 


1 target tvpe index V J 
p. pp ume period index ур тре P 
Parameter 

ТЕЙ length of time period p in davs Ур 
RECON, minimum number of days a target has to stay dead V J 
QTY, maximum number of targets } that can be reconstituted in 30 days V J 


Aggregated parameter 


J = index of the last time period considered for reconstitution. 


If a target of type j is Killed in time period PERUP,, or earlier, then there is sufficient 
ume available to reconstitute the target so that it once again will be available in period 
pti. Ihe parameters T/ME, and RECON, determine PERUP, according to the fol- 
lowing formula in HEAVY ATTACK: 


Let 

р+1 
00 1 if RECON, < >, TIME, ~ CEIL(0.5 x TIME) Vj,p<p<n 
p 0 otherwise p-p 


where the function CEIL rounds a real number to the next higher integer value. 


К; indicates whether targets killed in period p are eligible for reconstitution in 


period p and therefore: 
p 
PERUP,, — > Кр, م > مار و‎ 
ре! 
Note that always PERUP,, < р. 


Variables 


KILL, , number of targets tvpe } killed in time period p V 
REBUILD;, maximum number of targets of type j that are reconstituted 


as live targets in time period p+ 1 Үр я 
Conditions for Reconstitution 


A killed target of type j can be reconstituted if the following 4 conditions are true: 


l. at least a fraction of target j was destroyed in the previous or the current time pe- 
riod p, 


2. it has been dead for more than some defined time, 


3. the total number of targets being reconstituted has to be less than the total number 


үлүтртэн: 


of targets which exceeds the minimum dead time 


р РЕВСР., 
X REBUILD y < У КИШ, 
p=) p=) 


4. the maximum number of targets type |) which can be reconstituted at the end of 


VJ. po 


(A) 


(В) 


each ume period p s Si cit) 








| QTY, 
REBUILD TE qm x ГІМЕ, ү 
QTY, | 
where zg — Iepresents the reconstitution rate per day. 
This leads to the following submodel: 
max z= ) ? REBUILD;, 
ЦЭ, 
Bet: 
p РЕКЕР: 
ЕВО Е Vj p <r 
ru D 
| ӨТІ, | 
REBUILD; , < а^ Е De 


The interpretation of (A) is that the number of targets of type j rebuilt in period p or 


before cannot exceed the total number of targets that are killed during or before period 


PERUP,, . The interpretation of (B) is that the number of targets of type j rebuilt in 


period p cannot exceed a certain quantity depending on the length of period p and on 


the target tvpe. There are no targets reconstituted in the last time period p=n. 


ПІ. BOUNDS ON WEAPON CONSUMPTION 


A. INIRODUCTION OF A WEAPON CONSTRAINT 
A desired 1mprovement for the current HEAVY ATTACK model is to add an addi- 


tional constraint on the utilization of weapons inside the RAND - model. 


Two important facts should be recalled: 


1. For each sortie - target combination ( i, j ) and each weatherband there is at most 
one weapon which can be used. 


2. Averaging over all weatherbands is related to the probability that weatherband w 
might occur at the time sortie type i is assigned to target type j. 


Let the upper bound on weapon consumption be defined as: 


WP, total number of weapons of tvpe k available 


The required constraint for the consumption on weapons is then: 


v VN | | 
DS (Ура, х ВА Е Wan) 2 
7 и 


i 


Migere the Sum 1s over all ( 1, . w ) such that K 2 А, 


laj, W 


B. REVISED MODEL OF HEAVY ATTACK 

Reconstitution can be included in the RAND - model. Instead of considering re- 
constitution as a computational “bookkeeping” process, it can be part of the optimiza- 
tion. To accomplish this, it is necessary to define a new variable for the number of dead 
targets such that the time period as an additional dimension is represented bv a second 
subscript: 

D,, is the total number of targets of type j killed in time periods « p /ess the number 


of targets that are reconstituted during this time : 


--. 


pm 


Р. = 


JP (KILL; p ж REBUILD, 2) ү j. p 


л 


EN 


9, 


v 


V Jj 


vo 


УК 


The military worth of a target is also ume dependent: 
V. military worth of a tarzet type j in time period p 


Embellished Thesis Model (solved sequentially for p = 1, 2, 3,..., n) 


Max 2, = 2.05; x KILL,,) 


1 


SiG 


KILL, , = Пт Cj, S 27 Х рида X E. Ji 
i и 


where : f{... } is one of three functions discussed in the next chapter. 


) Љр 
p-i 
Dj, = 2 (KILL, , — REBUILD, „) 
р'=] 
9 PERUP, , 
REBUILD) < Пи 
pal р 


Эг”. 
j 


22 Б іх ده‎ x B, Ji < WP, 
[а 2 


where the sum is over all ( i, j, w) such that = 207 


0 € SX; 9 1, ل‎ 


RIEL чі 

2-0 ши 

О = REBUILD, , v 
where the upper bound on. REBUILD, 15 such that: 

QTY; | | 

221 йо = YJ 


The model was written in the General Algebraic Modeling System (GAMS, [Ref. 3). All 
optimization problems throughout the Thesis are solved with the nonlinear programing 
solver MINOS - Version 5.0 [Ref. 4. A database for 2 sortie-. 26 target- and 29 
weapon-tvpes was provided [Ref. 5] in order to compare the results by using three dif- 


ferent Objective functions, each over four time periods. 


IV. LINEAR VERSUS NONLINEAR MODEL 


In this chapter the derivation of the nonlinear objective function used bv RAND is 
given. In addition two alternatives are represented by introducing the Washburn- 
Equation and the linear case in which the number of killed targets 1s proportional to the 
number of assigned sorties. Each of the three objective functions 15 used in the model 
described in the previous chapter for sequentialy optimuzing sortie assignments over four 
time periods. In order to compare the effect of the three obedire UNCON я 


urement for the diversity of the allocated kill capability 1s defined. 


A. RAND EQUATION 
If A, represents the total number of killed targets of tvpe ) then the objective function 
used in the RAND - model can be derived from the differential equation: 


d K, K 
= | — ç x T (A) 





4 Х) j 


where А, = Ма хо паша тэг. 


i 


The differential equation (A) with the initial condition A(X, =0) — D, has the solution: 


Г, D, ; 
eu ED x [ — (1 — ох.) x E 


/ 
Instead of bounding K, by 
Does 
let AJLL, be the number of targets Killed in excess of D, : 
KILL, = K — D; 


so that 
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which leads to the final result: 


D 


. C; 
) — x x 
0) х (1 T ) 


KILL, = | 


B. LINEAR EQUATION 


A special case for the differential equation (A) appears when c, = 0: 





then 
d K; ET 
ЧА, 
which vields: 
A; = 4, Р, 
so that 
шэг 


or by using 
so that 


where the final solution represents the linear case: 


KILL, = X; 


illustrates the influence of Ше target parameter‏ لا 
с, оп the funcuon ATEL, — fU‏ 
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Figure 1. Influence of the target parameter c on the RAND-Equation: The sol- 
ution of the differential equation used in the RAND-model is graphically 
shown for 4 different target parameters c. 


The parameter c, has no direct physical motivation. The model considered in the next 


section also contains a single parameter. but the parameter can be motivated phvsically. 
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C. WASHBURN EQUATION 
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The Washburn - Equation [Ref. 6: p. 25] defines the differential E т in the fol- 
lowing way: di 
ак, 
от Probability ( attacking a live target ) 
д! 
or equivalently: 
dR { number of live targets } 
dA, {number of targets that look alive ) 
This leads to the differential equation: 
d K, Т, — К, 
j j j 
к » 
Ey б / D j 


where о, 1 а constant proportion of killed targets, which have the property to appear live 
to a potential attacker. 


The differential equation (B) with the initial condition K(X, — 0) 2 D, has the sol- 


Dj (I — x)x (K, — D) — 2, 
Kj = 1; x 15:12:55 сл Dm a 
y 


mune Ai instead of K such that: 


ution: 


KILL, — K, — D, 
leads to the implicit solution for the Washburn - Equation as: 


(1 — x)x KILL, — X, 
KILL, = (Tj = Dj) x (1 = e roD үр 
The difference between the two differential equations (A) and (B) for two different target 
parameters is shown in Figure 2 on page 24. Observe that for target parameter c close 


го 0 or 1 the Washburn-equation tends to behave similarly to the RAND-equation. 
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Target parameter o is denoted in the figure by с. 
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Figure 2. RAND- and Washburn-Diff. Equation with varied parameter с: The two 
differential equations are shown for 2 different target parameters c. Be- 
cause the solution of the Washburn-Equation can be given only in an 
implicit form, the differential equations are shown rather than their sol- 
utions. 


The influence of the three different objective functions on the RAND-model using 
the same input data is shown in Figure 3. 

The total worth of killed targets decreases with time for each objective function. 
The main reason for this is that in the first time period sorties are assigned to those 
target types for which the effectiveness is highest. When all targets are killed, sorties are 
then assigned in the following time periods to the remaining targets for which the effec- 
tiveness is less. As a result, more and more sorties need to be allocated in order to gain 
the same number of killed targets. The number of reconstituted targets available at the 
beginning of the second or third period is relatively small or even zero and can therefore 


be neglected at this point. Since the variation in the number of sorties and in the mag- 


nitude of the target values is too small to compensate for this effect, a declining trend 
in the objective function value over time for all three cases is observed. 

Note that the Washburn-Equation always yields a smaller value than the 
RAND-Equation. This follows from the fact that the Washburn-Equation declines faster 
than the RAND-Equation for the same target parameter c as shown in Figure 2. The 
linear equation is larger than either one. The most important difference is not in the 
absolute level of target value killed, but rather in the influence of the objective function 


on the distribution of sorties over targets. This subject is taken up in the next section. 
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Figure 3. Total Military Worth of Killed Targets: represented for each different 
objective function and each time period by the height of the respective 
block in the figure. 
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D. DIVERSITY OF KILLED TARGETS 

An important reason for USAF to use a nonlinear objective function is to avoid an 
undesired concentration of attacking sorties on a few targets. In analvsing the effect of 
the three different objective functions on the optimization, a measurement 15 needed in 
order to indicate how manv of the allocated sorties are spread over different targets. 


In information theory the function 


h(p) — Я (n x log = 


i 


where р (үд... Фу) and > = ] 
i 


is used to express the diversity or “entropy” of the probability distribution p = (p, j. 
Observe that A(p) = 0 when p concentrates all probability in one element. The maxi- 
mum possible value when p has n elements occurs when thev are all equal, in which case 
h(p) — logn. The diversity A(x) of an arbitary set {x,} of nonnegative members can 
be measured by simplv normalizing them so that they sum to | and then computing 


entropy: 


2,4) 
J 


The diversity of values A(x) gained from the same input data and model as used in the 
previous chapter 1s depicted in Figure 4. Since the number of targets n equals 26, the 


maximum diversity value will be 


Н(Х),д = 3.26 


Figure 4 makes it clear that the Linear objective function has a lower diversity value 


than the other two. This is to be expected, and in fact one of the main reasons for using 


a nonlinear objective in the first place was to avoid low diversity values. However, note 
that: 
]. The Linear diversity is not 0; that is, several target types are still attacked. 


2. None of the objective functions achieves complete (3.26) diversity. 


The differences emerge most strongly in period 3. Onlv 4 target types are attacked when 
the linear model 1s used, or 6 with the RAND-model. 16 different target types are at- 
tacked when the Washburn-equation is used; this is in keeping with the idea that the 
Washburn-equation is the most “non-linear” of the three (see Figure 2). The three 


models differ much less in period 1.2 or 4. 
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Figure 4. Diversity of killed targets for different objective functions: The height 
of each block illustrates to how many different target tvpes (out of 26) 
sorties are allocated at different time periods by using each of the three 
objective functions. 
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V. ALLOCATION OF SECONDARY WEAPONS 


A. COST-EFFICIENCY VERSUS KILL-EFFECTIVENESS 

Cost considerations are finished once SELECTOR has established the Preferred 
Weapon List. Although this list contains different tactics, ordered in terms of cost- 
efficiency, HEAVY ATTACK only uses the top one on the list which is feasible. The 
only time at which HEAVY ATTACK may proceed to a succeeding tactic appears. as 
mentioned before, when a weapon has been exhausted in earlier periods. 

As a second revision of HEAVY ATTACK, the model is changed to continue target 
attacks after the weapon type used by the most cost-effective tactic has been exhausted. 


using those weapons still on hand. 


B. A NONCONVEX CONSTRAINT 

The model discussed in the previous chapter requires that only the tactic on the top 
of SELECTOR's Preferred Weapon List can be used. Once the corresponding weapon 
type is depleted further attacks by that sortie type in that weatherband against that tar- 
get tvpe are impossible. The idea in this section 1s to relax this strict requirement to 
permit using whatever tactic is highest on SELECTOR 's list among those whose weapons 


have not been exhausted. 


Implementing this logic in the existing model requires a modification of the variable 
Эна 
SA 


TOU FW 


number of sorties of type i assigned to target of type J} which use 


tactic r in weatherband w 


The probability that all sorties of type i assigned to target of type j will attack the target 
in weatherband w has to be equal to the probability that weatherband w occurs at that 


time: 


r 


"s У у 
Л Хулм = PR, X (22:22 54р, u” 
7 2 


Upon these redefined variables for the number of assigned sorties, it 1s possible to de- 
termine the utilization of each weapon tvpe: 


let VEAP, be the consumption of all weapons of type k 


then ЕАР Е vk 
i j r w 


where the sum 1s over all ( 4, j, r, wj such e 


7 


In order to assign sorties using less cost-effective tactics, SX,,,,, must be 0 unless the 
weapon tvpes corresponding to all more cost-effective tactics are exhausted. The fol- 


lowing constraint will enforce this logic: 


r—] 


0 = SX, x 2, (FP, — ЕАР) У ії, тож (С) 


1 


r=] 


where К = A у. 


The above constraint requires that at least one of the two factors on the right hand side 
of the equation equals zero, so either no sorties are assigned (first factor zero) or else 
all more cost-effective weapons are exhausted (second factor zero). The constraint thus 
enforces the desired logic, but there is a disadvantage in using it. The disadvantage 15 
that the function on the right hand side of (C) is not only nonlinear (products of vari- 
ables are involved) but nonconvex. Without constraint convexity, there is no guarantee 
that the locally optimal solutions achieved by the '//NOS solver are globally optimal. 
There is some evidence, however, that globally optimal solutions are actually being at- 
tained. For one thing. employing constraint (C) always results in a higher objective 
function value than when only the most cost-efficient tactic 1s permitted. In addition, 
some experiments were performed where the improved model was changed into a linear 
model by linearizing the objective function at the optimal solution. The nonconvex 
constraint was then converted into a linear constraint by using integer variables. The 
optimal solution of this linearized model was identical to the solution gained by the 


nonlinear model with the nonconvex constraint. 
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C. REVISED MODEL 


Ihe mathematical model is solved sequentially for pa ыт 


Max z = 2 (M, x КИЛУ) 
7 


where Л, = 287 2 2211 


i и 
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where the sum is over all ( J, r, w} such that k = I S 
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where the upper bound on REBUILD, 


РЕШ 
RLBUILD, , | — 
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0-4 РАВ, < ИР, 


[he introduced relaxation. will be used in the further revision of HEAVY ATTACK 


رمه 
t.‏ 


considered in the next chapter. 


VI. GLOBAL VERSUS MYOPIC TIME OPTIMIZATION 


А. TIME-DEPENDENT MILITARY WORTH OF TARGETS 

When HEAVY ATTACK optimizes the allocation of sorties for each time period, 
it doesn't take advantage of the fact that the military worth of each target and each time 
period is Known prior to running the optimization. The decision, which target tvpe 
should be given a high priority to attack, 15 based on a comparison of military values 
of different target tvpes restricted to the current time period. Although military worth 
of a target is given as a function of time, HEAVY ATTACK doesn’t recognize the most 
favorable time for attacking a certain target type. This “myopic view” is caused by re- 
Stricting the optimization to the time interval covered by one period. 

It seems Worthwhile to consider an optimization covering all time periods at once. 
This “global” optinuzation is expected to spend resources even more effectively than 
before, so that the total sum of gained nuhtary worth of killed targets might become 
higher compared to sequential time optimization. In addition, it can be expected that the 
number and tvpe of killed targets in each time period will change. 

The third revision for HEAVY ATTACK as presented in this chapter doesn't require 
major changes to the previously discussed model. A subscript for time is added to the 
Able Sew: 
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isl WP 


number of sorties of type i assigned to target type j by using tactic 


tvpe r in weatherband w and in time period p 


The resources on sorties available needs to be defined as a function of sortie tvpe and 
time: 
S maximum number of sorties type 1 available in period p 


PROP,, proportion of S,, that can be assigned 


Computing time increases with the number of time periods covered. 


B. GLOBAL MODEL 
The mathematical model is shown below. The realization of this model in GAMS, 


including all inputs, 15 given in the Appendix. 
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Where the upper bound on REBUILD,, is such that: 
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C. RESULTS AND COMPARISONS 

Ihe above model was too large to be run in GAMS on available computer equip- 
ment at reasonable cost with the same size of input data used previouslv. Therefore the 
number of target tvpes were reduced from 26 to 13. Other efforts were also made to 
decrease required computing time. | 

Table 10, Table 11 and Figure 5 compare the results of the global and mvopic se- 
quential optimizations. The global optimization achieves more target value killed; the 
percentage gain for the global approach 15 (1358.0 - 1123.0) 1123.0 - 20.9 %. Com- 
paring the target values of target tvpe 5 and 27 over all 4 periods shows that the highest 
target value occurs in period 3. The global optimization realizes this fact by destroving 
all available targets at that time. While both target types, especially target tvpe 5. have 
a relatively high target value in the first time period, most of these targets are therefore 


killed by mvopic optimization in the first period. 


Table 10. NUMBER OF KILLED TARGETS: The table shows the number of 
killed targets achieved by sequential and global optimization as well as 
the ср target value for each time period. 
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Table 11. MILITARY WORTH OF KILLED TARGETS: gained bv sequential 
and bv global optimization is given for each time period and as a total 
sum. 


Optimization Optimization 


Total Worth 1 
of Killed Targets 1123.0 1358.0 
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WORTH OF KILLED TARGETS 


0 


Figure 5. Distribution of Military Worth of Killed Targets: The height of each 
block represents the numerical value given in Table 11 depending on the 
time period and on the kind of optimization used. 
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Both the global and the mvopic models utilize secondary weapons. Figure 6 shows 
weapon usage in the global model. Note that weapon tvpe WP7 is used extensivelv in 


situations where more cost-effective weapons are exhausted. 
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Figure 6. Allocation of Secondary Weapons: The height of each block represents 
the number of weapons utilized by the global optimization. A significant 
number of weapon tvpe WP7 18 used bv tactics Gf order r=3. 1015 15 
only possible when weapons used by tactics of order r=1 and r=2 are 
exhausted. 


A more detailed report of the solution is given in the SOLVE SUMMARY of GAS in 
the Appendix. 
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VIE- CONCLUSIONS 


In the first revision of the current HEAVY ATTACK model, a weapon constraint 
is added and three different objective functions are compared. The objective function 


best used in the model depends on the priorities of the user: 


1. Using a linear objective function instead of a nonlinear one has the advantage of 
simplicity and consequent computational efficiency. A. disadvantage is a less dis- 
persed allocation of sorties to different targets. 


2. Using the Washburn - Equation instead of the RAND - Equation has the advan- 
tage of using a well defined target parameter. The dispersion of attacked target 
types might be somewhat less influenced due to changes in the input data. 

In the second revision the current philosophy of using the most cost-efficient tactic 
is relaxed such that less cost-efficient tactics can be utilized within a time period. With 
this revision, tactics not at the top of the Preferred Weapon List (SELECTOR output) 
can be utilized if all more cost-effective tactics are infeasible due to Weapon exhaustion. 
This revision is particularly important when there 1s a small number of time periods. 
since the same capability already exists between time periods. 

The third revision replaces sequential optimization (current practice) with global 
optimization. The comparison between sequential and global optimization by using the 
same input data shows a qualitative difference in the achieved results. There 1s a definite 
indication that sequential time optimization tends to achieve military. success in the be- 
einning of the war bv sacrificing the potential for later success. Global optimization 
tends to husband weapons and even targets (in cases where target value increases with 
ume) for later periods in the war. An argument for global optimization can be based on 
the fact that it is more efficient in killing targets with large military values. On the other 
hand, it could also be argued that sequential optimization is more likely to imitate what 
Will actually happen, “optimal” or not. In any case, if global optimization is used, then 
the distribution of the value of destroyed targets seems to be much more time dependent 
than is recognized by the current method of sequential optimization. 

All revisions introduced in this Thesis result in gaining of more military. worth. 
USAF's general objective is to determine their future need of weapons rather than to 
maximize the military worth of killed targets. With the revisions described above, utili- 


zation of weapons plays a more important and direct role in the optimization, especially 
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when more than one tactic is considered. The developed models are intended to provide 
the necessary structure to embellish HEAVY ATTACK for this purpose. 


40 


APPENDIX GLOBAL OPTIMIZATION MODEL 
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is allowed to be reconstituted 


QTY( J) 


LOOP((J,P), 


BERUBBSUP) — SUM(BPSCORDCPP) БЕ ORDCPO)OSISORBEEONCJ) LT (SUNK PPPS 


(I O TOTI PP) DLESCORDCP)T10).AND (ORDCPPBEOUERMSRESPP)) ),Т1МЕ(РРР)) 


E 


= СЕТИО. 5 2 Өө »))) 


ус 


v9 ev (v 0C 2V 4V dV м كل كي كل كل‎ ЭЭ كل ول كن 6ل‎ Hv dv كل‎ UV V ZV dv dV QV كي‎ v dV VV VV Vv 4A dv de VV 4e dA dV de dV VV VV Vv dV VV d dV d v 


Begin of aggregated INPUT DATA Үе 


Эс 
Ус 


sc 


fo nl o ufo‏ اوو اہ ТУ іы Усы‏ لن ی ان یگن ی لن وان یگن ی کن ےن یکن یکن ی ن مان م کن وگن مان كما أي أ وكا وك وإ ىأ مإ ملع أ وأ ولو ولك وك وليك وكا لوأ مأو لو ملو مأك 
е еъ съ Фу Фу 48 4S 4S S 4S0 ER PR EL GD ER FR AY EL FH TV ER FB 4х EH Фъ 04S 4S 6% 4% 4% 4% 04S 4S Фу 4% 6% 4S 4S IS съ ду чу чу 45 4 eves as 49 #١ » 27% 4»‏ 
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Number of Targets type J killed by one Sortie type I 


OD1.WB2 011.183 OD1.WB4 0D1.WB5 +6 


دم دم око КМ СК КО FS S‏ كن كن كن كن ذم دم دم £O‏ 


004 


‚ОШО 
327 
“2716 
. 034 
. 044 
2209 
21297 
1905 
„216 
‚042 
. 741 
ELI 
. 428 
1282 
. 218 
26 
. 034 
. 044 
2209 
39065 
59250 
135 
. 042 
2212 
2092 
1211 


OD3 


2195 
2299 
. 076 
2077 
2220 
«205 
. 284 
22225 
1203 
2067 
20/192 
2092 
2150) 
290 
7065 
.076 
2000 
2020 
„260 
509 
. 468 
25 
7067 
2215 
21:33 
20908 


OD2 


1156 
. 418 
23:20 
2092 
. 010 
„269 
2222 
. 468 
252 
-072 
: 
212) 
зи 
. 241 
205 
5120 
1092 
2010 
„216 
(52 
2383 
206 
2072 
29139 
.064 
1967 


OD1 
2159 
1305 
2063 
Eoo 
028 
2216 
2056 
. 343 
2273 
. 134 
4933 
2138 
226 
. 7 
2702 
-083 
-061 
. 028 
«42955 
2955 
251 
. 205 
. 134 
10572 
А157 
52 


Weaponload Array for each set <i j r w? 


су FS ÉS GN ON O. ON G. جم‎ го SS ON دحم دم دمح دمح دم دح‎ 


Гг كح ذنم دم جم كل‎ ON O. O. O. IS го ко СУ MN دم دمح‎ $S 


£ ON O МЮ МЮ О О G, ON GO. S O tS @ FS оо го го N 


СО ON O IÑ IS O. Os ON O O ES O FS ƏO ES O ES NO ES N 
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OD1. WB1 


ооо осо ооо 


TGs 

TGE 

TG10 
ТБ МІ 
1612 
TGS 
TG14 
TG22 
TG24 
T6025 
o2 
1629 
TG34 
TG5 

TG8 

TG10 
TG11 
TG12 
TG13 
TG14 
TG22 
TG24 
TG25 
T6025 
TG 29 
TG34 


TABLE Ес 


АС1. 
AC1. 
АС1. 
АС1. 
АСТ. 
АСТ. 
АСТ. 
АС1. 
АС1. 
АСІ. 
АС1. 
АС1. 
АС1. 
АС2. 
AC2. 
AC2. 
AC2. 
AC2. 
АС2. 
АС2. 
АС2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 


TABLE  B(I,J,R,W) 


АС1. 765 
АС1. 768 


ШОО 
TO TE 
SII 
/45413 
. 4 
I622 
. TG24 
71525 
‚ [627 
21629 
. 1634 
2» 

IGG 

LODS 
С 
Б 
21019 
. TG14 


с г гг г ос" 


002. КВ6 


Хосоозозосоаоосососососеасссссо 


003. КВ6 


نت دخ دخ وح روح يح يم يم جم ب نب ذم نب بن ب تم دم ديم ات 


NON AAD 


OD2.WB2 OD2.WB3 OD2.WB4  OD2.WB5 


OV ON O O O NS O S O O O S OA N O O O O. O O O O O O O N 


OD3.WB2 ODS. WBS OD3.WB4 OD3.WB5 


O O FS IS tS O O O kt. O O kb. PS S tS SS tS: tS ES O 


ОУ О КЮ O. O. 


ON ON ON O O ífS O NS O O O sS OA N O O O O. O O O O O O O FN 


O ON NN OOON OON NN NNN N N O 


£- ON O FS O. O. 


со O O. O O N O S O O O iS O O O O O O. O O O O O O O N 


OONNNOON NOON NNN NN N N O 


ооо ON O. 


C2 CO CO CO (e) т о ОЭ оо оО Ооо о © О ӘӘ 


OOONNOON OT OONNN ON N N O 
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СО СОО СО О со 


OD2.WB1 


SO GO Фо СО СО СО СИК) СО © O СЗ С) СО КО СО Со Со ко со кое تك‎ 


орз. ВІ 


сос си со со CO O О СО 3 3: 5ك‎ 5 CO © © 


1622 
TG24 
1222) 
1027 
Ш529 
TG34 


AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 


AC1. 
AC1. 
AC1. 
AC1. 
АС1. 
АС1. 
АС1. 
АС1. 
АС1. 
АС1. 
АС1. 
АС1. 
АС1. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
ACA 
AC2. 
AC2. 


233 
334 
695 
336 
557 
299 
209 
340 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 
251 
252 
253 
354 
5255 
5256 
257 
358 
Dog 


OM, ©‏ تك © تك 


OD4.WB2  OD4.WB3 0104.МВ4 OD4&.WB5 014.6 


ОУ КОСО СО СО OCOD Со С ON OGO O СО ك2 5 لك‎ OO © GO со КОСО 


оооомо 
сос O O O i O 


O S O O O O O. O O O FS O ƏO O O bk O ƏO QO. O OQ. O O kb. G. O. 
O S O O O O ОӘО O O kS O O O O kS O ƏO GO. O QO. O O bF. O. O. 


* For each sortie-target combination the weapon type K of order R 


г 


O 
о ےہ یہ یہ‎ с о о о (ліло о г О 


= 


М 


ON O O O O O 


ош”юшссососоуүбоссоус сс сэ. сс бос Ээ со со Өс 


сос‏ حك 


ооо оо о во о ә оооло: 


004. WB1 


ооо 


СО ОО Со Оо © فكت كت‎ со» сс هك‎ ээ CO CO СО СО СО СО Со CEOS SE © 


22 
TG24 
1525 
1827 
216) 255 
TG34 


if it is possible to use this weapon 


OD3 


л хм Оз ФУ ЧА SI л <J <J >J мы м м 


ы 


O 
Мә 


(л Co Co FP Co Ui Co Ui Ui Ui Ui Ui Os (CJ 


قم 
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— 


G U OU т Ui Ui Co Ui Ui Ui Ui Ui Un С 


© 


ho 


AC2. 
AC2. 
AC2. 
AC2. 
AC2. 
AC2. 


AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC]. 
АС1. 
АС1. 
АС1. 
АСЫ 
АС1. 
АС2. 
АС2. 
АС2. 
АС2. 
АС2. 
АС2. 
АС2. 


АС2. 
АС2. 
АС2. 
АС2. 
АС2. 


* is given 


155 

TG8 

TG10 
CII 
15012 
153 
TG14 
1022 
TG24 
ШЕ/2 5 
11527 
ШБ 2 9 
TG34 


AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
AC1. 
АСІ. 
АС1. 


TABLE, КЕШЕГІ, Л2К) 


J 


229 
320 
SOT 
292 
295 
394 
295 
226 
397 
295 
299 
400 
401 
402 
403 
404 
405 
406 
407 
408 
409 
410 
411 
412 
413 
414 
415 
416 
417 
418 
419 
420 
421 
422 
423 
424 
425 
426 
427 
428 
429 


431 
432 
433 
434 
435 
436 
437 
438 
439 
440 
44] 


443 
444 


ж 


АС2. 765 2 1 2 1 
AC2. TG8 1 5 1 5 
AC2. TG10 5 S 7 18 
AC2. TG11 5 5 7 5 
АС2. ТС12 5 5 7 5 
AC2. TG13 1 5 1 5 
20224 12 12 12 1 
АС2. 7622 1 5 5 1 
AC2. TG24 1 1 5 1 
AC2. TG25 24 24 24 24 
AC2. TG27 1 3 9 3 
AC2. TG29 3 1 3 7 
AC2. TG34 12 1 1 1 
еее 
End of INPUT DATA * 


سا 
* 


“... 


ae 
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Definition of Sortie Variable 


SX(I,J,R,W,P) describes the number of sorties type I assigned 
to a target of type J carrying any weapon feasible for tactic R 
and weatherband W and in time period P 


SACI TIRS W, R) 


POSITIVE VARIABLES 


Initial Values for Variables 


SX.L(I,J,R,W,P) = 0 


Declaration of variable EXPO(J,P) 


POSITIVE VARIABLE EXPO(J,P) 


Declaration of Kill Variable 


POSITIVE VARIABLE KILL(J,P) 


Declaration of Variable D(J,P) 


POSITIVE VARIABLE D(J,P) 


Declaration of Variable for cumulative weapon consumption 


POSITIVE VARIABLE WEAP(K) 


Upper bound for variable Weapon Consumption 


= WP(K) 


49 


WEAP. UP(K) 


ть ом 


зоог 


500 


501 ** Declaration of variable for number of targets been reconstituted 
502 

503 POSITIVE VARIABLE REBUILD E) 
504 

505 ** Upper bound for variable REBUILD 

506 

507 REBUILD.UP(J,P) = QTY( J) ü TIME 2. : 
508 

509 

510 

511 ** Variable definition for objective function 

312 

513 VARIABLE Z ; 
514 

515 

516 EQUATIONS 

217 

518 KILLVAL maximize the value of destroyed targets 

519 KILLNL(J,P) determines the number of killed targets 

529 EXPONENT(J , P) evaluates the values of the exponential terms 
521 DEADTG J P) determines the number of dead targets 

522 KILLCON(J,P) constraint the number of killed targets 

228 REGCOCOD TP) constraint the max. number of targets for reconst. 
524 SORTCON(I,P) constraint the number of allocated sorties 

525 WEAPCONSUM(K) determines the consumption of each weapon type 
526 SELECT(I,J,R,W) decides if next weapon on list can be used 

927 ЕЕ. р» ensures that all weatherbands are covered prop. ; 


ЗРО ТОА 


5 Z =E= SUM((J,P);V J p) КТ УР) ; 


среди Керъл TT 


529 KILL(J,P) =E= ( (T(J)/C(3J)) = тт р. D i 


540 EXPONENT(J,P).. 
542 EXPO(J,P) sE- EXP( ((-C(J))/T(J)) * SUM((I,R,W)SB(I,J,R,W), 


БАЙ E(I,J,R) * SX(I,T,;R,W,P)SB(I,J EDD 


547 DEADTG(J,P).. 


549 D(J,P) =E= SUM(PPS(ORD(PP) LT ORD(P));,KILL(J Pp) = REBUILT ee DER 


S52, KILLCO J p m 
554 KILL(J,P) =L= T(J) =SP(J,P) : 


556 





RECCON(J,P).. 


SUM(PPS(ORD(PP) LE ORD(P)),REBUILD(J,PP)) =L= 


SUMNENESSORDCPP) LE PERUP(J,P)),KILL(CJ,PP) ) ; 


ЗӨ КТСӨМСТ,Р).. 


SUMCCJ,R,W)$B(I,J,R,W), 


ВВМ) =I PROP(I,P)SSESCI,.P) ; 


WEAPCONSUM(K).. 


WEAPCK) -Е- SUMQC(I,JQRSW,P)SC CORB(K) EQ WPTYPE(I,J,R)) AND 


(B(I,J,R,W) NE 0) ),B(I,J,R,W) * SXCI,J,R,W,P)$B( I,J,R,V)) ; 


SELECT( I,J,R,W)SB(I,J,R,W).. 


О = Е SUMCPSSRCEIIJSR,W,POISBCI,J,R,W)) * 
SUMCCK,RP)S( CORDCRP) LT ORD(R)) AND 


В hp W) NE 0) ЖИН СОКБОК) ЕО WPTYPE(IjTSSRP)) ), 


5 


5 


(WP(K) - WEAP(K)) ) 


values 


DISTR(I,J,W,P)SSUM(R,B(I,J,R,W)).. 


SUM RAS T.J R W P) 5B I, ,J,R.,W)) -Ез РКСУ) 7 


SUNCCR ВВ В. Е, РР), ХЕ, P b P P)SB( (TI, ,J,R,WPP)) ; 


MODEL AIRATTACK /ALL/ 


ve 


Limit for number of iterations 


OPTION ITERLIM = 1000 , LIMCOL = 0 , LIMROW = 0 


OPTVONVS@EPRINT — ORPOSSSTSUSIT- OFF 


SOLVE AIRATTACK USING NLP MAXIMIZING Z 


T The following statements represent the solution 


PARAMETERS 


557 
558 
559 
560 
56l 
562 
563 
564 


. 565 


566 
567 
568 
509 
SM 
S 
572 
573 
57& 
575 
3/6 
577 
58 
579 
580 
581 
307 
303 
584 


606 
607 
608 
609 
610 
611 
22 


613 





614 KILLTG JTE) number of targets J killed in period P 

615 OBJECTIVE(P) Objective Function Value | 
616 KILLEOT JITE) potential Kill-Capability (target-type vs period) 
617 OPSORTIE(I,J,R,P,W) number of optimal sorties 

618 SORTIE C J P D) number of sorties I assigned to target J in period P 
619 200281 ل‎ number of weapons (sortie , target and weapon type) 
620 WPCONS(R,K) number of weapons (tactic vs weapon-type) | 
621 WEAPON( J,K) number of weapons (target vs weapon-type) : 
622 

65 

624 KILLETG(J;, P) = KILCLAECCJS EP) ; 
625 

626 ОВЈЕСТІУЕ(Р) = SUMCJ,V(J,P) * KIL БОЕ : 
627 

628 КТ РОТС3 ,Р)у“52э201С0ЕК2 (5 Е ЕГЕ 

629 

630 БС1,4,Ку ме 5Х: нео У 0? : 
631 

632 WEAPONCJ,K) = SUMCCI,R,W,P)SCORDCK SEO Pie Gk ae 

25% 

634 B(I,J,R,W) = SX. L(I, J. RIV P eee 
635 

636 WPCONSC(R,K) > 5 PST 

637 

638 (ORD(K) EQ WPTYPE TI J.R ANDA RO TS TONS EUM 

26 

640 B(I,;J,R,W) ~ SX. LCI; J, R 7H " 
641 

642 ОРБОКТТЕСТЭ2 КР (252522 1219 0 ; 
643 

644 SORIIEC J P T) = SUM((R,W)RSX. L(I,J,R, ,Wp)) ; 
6&5 

6&6 

6&7 

646 OPTION OBJISGTDVE:2 DISPLAY OBJECTIVE 


DISPLAY KILLTG 
DISPLAY KILLPOT 
DISPLAY OPSORTIE 


649 OPTION KILLTG: 1: 1:1 ; 
DISPLAY SORTIE U 
$ 


СЭ ОРТТО КЕ 2-1 
5517 0211010 0 0111: 1: 
ООО в ОКЕ ШЕ: 1: Ще 2 
Ша DISPLAY WPCONS 
ТІ DISPLAY WEAPON 
DISPLAY WEAP.L 


DISPLAY REBUILD.L 


655 "Or PPO I PWP SON: 1: 

654 OPTION WEAPON: 1: 

бээ? OPTIONS WEAR? 1 

63C OPTION REBUTED: SIM 


ч э we we We we we фе фә NA 


COMPILATION TIME = 2.140 SECONDS 


MODEL STATISTICS SOLVE AIRATTACK USING NLP FROM LINE 607 
MODEL STATISTICS x 
| 


| 


| 


BLOCKS OF EQUATIONS 10 SINGLE EQUATIONS 932 


BLOCKS OF VARIABLES 7 SINGLE VARIABLES 1289 
NON ZERO ELEMENTS 2758 NON LINEAR N-Z 1889 
DERIVATIVE POOL 3 CONSTANT POOL 61 
CODE LENGTH 15943 

GENERATION TIME = 82-0 SECONDS 


EXECUTION TIME 67.680 SECONDS 


SOLUTION REPORT SOLVE AIRATTACK USING NLP FROM LINE 607 
SOTE SUMMARY 
MODEL АТКАТТАСК OBJECTIVE 2 
ШУРЕ МІР DIRECTION MAXIMIZE 
SOLVER MINOS5 FROM LINE 607 
pest SOBVER STATUS 1 NORMAL COMPLETION 
weeks MODEL STATUS 2 LOCALLY OPTIMAL 
were OBJECTIVE VALUE 55. 0172 
RESOURCE USAGE, LIMIT 64. 179 1000. 000 
ITERATION COUNT, LIMIT 639 1000 
EVALUATION ERRORS 0 0 
NINOS “VERSION 5.0 APR 1984 


COURTESY OF B. A. MURTAGH AND M. A. SAUNDERS, 
DEPARTMENT OF OPERATIONS RESEARCH, 
STANFORD UNIVERSITY, 
STANFORD CALIFORNIA 94305 ОП. Вай, 


WORKNSPACE NEEDED (ESTIMATE) -- 104191 WORDS. 

WORK SPACE AVAILABLE -- 134740 WORDS. 

(MAXIMUM OBTAINABLE -- 288878 WORDS. ) 

гэг ОРГТЇТМАГП SOLUTION FOUND 

MAJOR ITERATIONS 22 

NOR kG NORM РІ 5.752Е-06 

TOTAL USED 6524 UNITS 

BNNOSD TIME 56 7 (INTERPRETER: DA) 
c REPORT SUMMARY : 0 NONOPT 


0 INFEASIBLE 
O UNBOUNDED 
0 ERRORS 


648 PARAMETER OBJECTIVE OBJECTIVE FUNCTION VALUE 


БЕК 251959 БЕК2 333784 PER3 674.04, РЕК4 98.84 


649 PARAMETER KILLTG NUMBER OF TARGETS J KILLED IN PERIOD P 


PERI PER2 PER3 РЕКА 
По 0.5 i 18.0 2.0 
TG8 13.0 
ТОПО 26.3 
Тт 9.6 
TG13 272 1 с; 0.1 
Teale. 2.0 
T? 6.0 
тера 1.5 1.3 
TOP 26.6 1-1 
'T627 2170 
TG34 18.0 
---- 650 PARAMETER KILLPOT POTENTIAL KILL-CAPABILITY (TARGET-TYPE VS 

PERIOD) 

PER1 PER2 PER3 PERA 
TGS 0.5 lea 20.1 2.5 
WES 113207 
TG10 29.0 
mST 9.8 
ТОЗ 2.4 277: (NO 
ТЕ PT 
T. 6.7 
TG24 1.9 3.3 
1625 2879 1.3 
210727 36.1 
TG34 20.1 


651 PARAMETER OPSORTIE NUMBER OF OPTIMAL SORTIES 


INDEX 1 = AC1 INDEX 2 = TG8 
WB2 WB3 МВА WB5 WB6 
0D3. PERI 0222 1 0.6 06 520 


INDEX 1 = AC1 INDEX 2 = TG10 


WB2 WB3 WB4 WB5 WB6 
OD1. PER4 4.2 29.4 14.7 14.7 9273 
003. РЕК4 47. 1 





INDEX 1 - AC1 INDEX 2 = TG11 


54 


WB6 


WB6 


WB6 


2 


WB6 


109.4 


WB6 


WB5 


WB5 


WB5 


Сэл СО 
м мо м 


WB5 


WB5 


WB4 


WB4 


O N 
“IO 


WB4 


2:09 


WB4 


WB4 


сл 


ел 


WB3 
14.0 
То 
WB3 
1.4 
T9795 
WB3 
14.0 
1-3 
TG27 
WB3 
5.4 
TG34 
WB3 
55 
TG5 
WB3 
05 
SO 
11.4 
1.4 
TGS 
WB3 
60 
[610 
WB3 
2159 
Та | 
WB3 


WB2 
20 
INDEX 2 
WB2 
0-2 
INDEX 2 
WB2 
20 
0. 2 
INDEX 2 


WB2 


0.5 
INDEX 2 

WB2 

0.9 
INDEX 2 


WB2 


АС1 


АС1 


АС1 


АС1 


АС2 


qa 3-2 
925020522 


16 
02 


INDEX 2 
WB2 
019 

INDEX 2 
WB2 
За 

INDEX 2 


WB2 


AC2 


AC2 


AC2 


ODIT PERI 


INDEX 1 = 


oil PERS 


INDEX 1 = 


CPL PERL 
DUE PERS 


INDEX 1 = 


ЕМЕ PERS 
ODS TERS 


INDEX 1 = 


ODI PER3 


INDEX 1 = 


ODIT PERI 
LE PER2 
OPT PERS 
001. РЕКА 


INDEX 1 = 


ОРЕК 


INDEX 1 = 


OD3. PER4 


INDEX 1 = 


ODI EERI 0.4 21) >, 1 14.6 


INDEX 1 = AC2 INDEX 2 - TGYIS 

WB2 WB3 ИВА WB5 WB6 
OD Perak OFZ 77 0.8 0-9 8.4 
0р1. РЕКА 1 31501092 ШЕСІ 6.21 - 2 CONSE 056 
INDEX 1 = AC2 INDEX 2 - TG14 

WB2 WB3 ИВА WB5 WB6 
ОВ ВЕК 0.4 072 222 212 
OD STERI 0202 
INDEX 1 = AC2 INDEX 2 = TG22 

WB2 WB3 WB4 WB5 WB6 
0р1. РЕКА 0.5 S 1 9 9 Lo 
INDEX 1 = AC2 INDEX 2 - TG24 

WB2 WB3 WB4 WB5 WB6 
ШЕШЕРЕЛКЗ 02 10 
0р1. РЕКА 9.3 2-43 № 1 1.1 Ll 
05 ВЕК 05 ل‎ 221 
INDEX 1 = AC2 INDEX 2 = TG25 

WB2 WB3 WB4 WB5 WB6 
ODE ER 2 2 LO? 2-1 5 1 51. 
INDEX 1 = AC2 INDEX 2 = TG34 

WB2 WB3 WB4 WB5 WB6 
DI PERS 233 ИЫ 1::1 11.4 
(DESDE з 0. 3 
s 652 PARAMETER SORTIE NUMBER OF SORTIES I ASSIGNED TO TARGET J 

INSCEERIODEE 


PER1. AC1 PER1. AC2 PERZ. ACI PER2. AC? PER3. AC1 PER3. AC2 
ше, 2 4.3 01:7 
TG8 9.8 42. 
TG11 22) 20: 
T6179 12: 
ТС 14 Эг 
TG24 7.4 
TG25 10079 11521 223 


LORS 


\O © М‏ م مم 





828 
46.5 15-23 
PER4. AC1 РЕКА. АС2 
TOS 
2100 15623 
0779 
26-1 
[o 
653 PARAMETER WPCONS NUMBER OF WEAPONS (TACTIC VS WEAPON-TYPE) 
МР1 МР2 КРЗ WP5 WP7 TIZ 
959820 550628 53079 556 25.2 
022 
20 4.7 12.4 423.5 
WP46 
450.0 
654 PARAMETER WEAPON NUMBER OF WEAPONS (TARGET VS WEAPON-TYPE) 
WP1 WEZ WIES WP5 WP7 WP12 
586. 8 
256.6 16.7 
3258 406. 9 
241. 2 
1 Z 20-7 
12.4 
160.0 
OS 12.4 
29820 
20 271420 64.1 
WP46 
450.0 


655 VARIABLE WEAP. L 
600.0, WP2 586.8, WP3 511.9, WPS 600.0, WP7 423.5 
76.4, WP46 450.0 
656 VARIABLE REBUILD. L 
PER1 PER2 PER3 


9.5 1221 2.0 


1527 
TG34 
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БУО 
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1120272 
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OD1 
OD2 
OD3 


OD1 


Т65 
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TG10 
ШӘЛІ 
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TG14 
1622 
1624 
1627 
TG34 
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МР 
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165 


ШӘДІ 
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ло 
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EXECUTION TIME 22.580 SECONDS 
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